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Abstract. We give a partly new proof of the fluctuation bounds for the second class 
particle and current in the stationary asymmetric simple exclusion process. One novelty 
is a coupling that preserves the ordering of second class particles in two systems that are 
themselves ordered coordinatewise. 



1. Introduction 

The asymmetric simple exclusion process (ASEP) is a Markov process that describes 
the motion of particles on the one-dimensional integer lattice Z. Each particle executes a 
continuous-time nearest-neighbor random walk on Z with jump rate p to the right and q to 
the left. Particles interact through the exclusion rule which means that at most one particle 
is allowed at each site. Any attempt to jump onto an already occupied site is ignored. The 
asymmetric case is p ^ q. We assume < q < p < 1 and also p + q = 1. 

In this paper we give a partially new proof of the fluctuation bounds for ASEP first 
proved in [TJ and [2]. For a more thorough explanation of the context and related work we 
refer the reader to these earlier papers. The present paper has two novelties. 

(i) Both the earlier and the present approach are based on identities that link together 
the current, space-time covariances, and the second class particle. In [2] we proved these 
identities with martingale techniques and generator computations. Section [2] of the present 
paper shows that these identities are consequences of not much more than simple counting of 
particles. These arguments are elementary and should work very generally. Without much 
effort we extend the identities to more general exclusion dynamics that allow bounded jumps 
and rates that depend on the local configuration. However, product invariant distributions 
remains a key assumption. 

(ii) In pQ we used sophisticated couplings introduced in [I] and some complicated esti- 
mation to bound the positions of certain second class particles. In Section [3] of the present 
paper we introduce a coupling for ASEP that keeps the second class particle of a denser sys- 
tem behind the second class particle of a comparison system. Since the macroscopic speed 
of a second class particle is H'(p), the slope of the flux at the given density p, concavity of 
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the flux suggests that second class particles travel slower in a denser system. Since this is 
what the coupling achieves, we think of this coupling as a form of microscopic concavity. 

Fix two parameters < q < p < 1 such that p + q = 1. We run quickly through the 
fundamentals of (p, g)-ASEP. We refer the reader to standard references for further details 

mm- 

Definition and graphical construction. ASEP represents the motion of particles 
on the integer lattice Z. The particles are subject to an exclusion rule which means that 
each site of Z can contain at most one particle. The state of the system at time t is a 
configuration n(t) = (r]i(t))i & i G {0, 1} Z of zeroes and ones. The value rji(t) = 1 means 
that site % is occupied by a particle at time t, while the value rji(t) = means that site i is 
vacant at time t. 

The motion of the particles is controlled by independent Poisson processes (Poisson 
clocks) |JV i_>i+1 ,iV i->i_1 : i £ Z}. Each Poisson clock jV w * +1 has rate p and each N 1 ^ 1 ^ 1 
has rate q. If t is a jump time for N l ~ n+1 and if (77^ (t — ), ?7i+i (t — )) = (1,0) then at time t 
the particle from site % moves to site i + 1 and the new values are (rji(t),rji + i(t)) = (0, 1). 
Similarly if t is a jump time for iV 1 ^ 1-1 a particle is moved from i to i— 1 at time t, provided 
the configuration at time t— permits this move. If the jump prompted by a Poisson clock 
is not permitted by the state of the system, this jump attempt is simply ignored and the 
particles resume waiting for the next prompt coming from the Poisson clocks. 

This construction of the process is known as the graphical construction or the Harris 
construction. When the initial state is a fixed configuration 77, P v denotes the distribution 
of the process. 

We write n, u>, etc for elements of the state space {0, 1} Z , but also for the entire process 
so that 77-process stands for {%(£) : i £ Z,0 < t < 00}. The configuration S{ is the state 
that has a single particle at position i but otherwise the lattice is vacant. 

Invariant distributions. A basic fact is that i.i.d. Bernoulli distributions {^ p }pe[o,ii 
are extremal invariant distributions for ASEP. For each density value p € [0, 1] , v p is the 
probability measure on {0, 1} Z under which the occupation variables {77^} are i.i.d. with 
common mean J rji dv p = p. When the process rj is stationary with time-marginal u p , we 
write P p for the probability distribution of the entire process. The stationary density-p 
process means the ASEP 77 that is stationary in time and has marginal distribution rj(t) ~ 
u p . 

Basic coupling and second class particles. The basic coupling of two exclusion pro- 
cesses 77 and u means that they obey a common set of Poisson clocks {JV*-** +1 , JV* - ** -1 }. 
Suppose the two processes 77 and i] + satisfy t/ + (0) = 7/(0) + <Wq) a t time zero, for some 
position (5(0) € Z. This means that r]f(0) = ^(0) for all i 7^ Q(0), Vq(o){Q) = 1 an d 
7 ?Q(o)(0) = 0. Then throughout the evolution in the basic coupling there is a single dis- 
crepancy between n(t) and r] + (t) at some position Q(t): rj + {t) = 7](t) + 8qm. From the 
perspective of n(t), Q(t) is called a second class particle. By the same token, from the 
perspective of rj + (t) Q(t) is a second class antiparticle. In particular, we shall call the pair 
(77, Q) a (p, q)-ASEP with a second class particle. 

We write a boldface P for the probability measure when more than one process are 
coupled together. In particular, P p represents the situation where the initial occupation 
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variables 77^ (0) = ^ + (0) are i.i.d. mean-p Bernoulli for i 7^ 0, and the second class particle 
Q starts at Q{0) = 0. 

More generally, if two processes 77 and to are in basic coupling and to(0) > f?(0) (by which 
we mean coordinatewise ordering tOi(0) > 77^ (0) for all i) then the ordering to(t) > n(t) holds 
for all < t < 00. The effect of the basic coupling is to give priority to the 77 particles over 
the uj — r\ particles. Consequently we can think of the ^-process as consisting of first class 
particles (the r\ particles) and second class particles (the to — 77 particles) . 

Current. For x £ Z and t > 0, J x (t) stands for the net left-to-right particle current 
across the straight-line space-time path from (1/2,0) to (x + l/2,t). More precisely, J x {t) = 
Jx{t) + — J x (t)~ where J x (t) + is the number of particles that lie in (— 00, 0] at time but 
lie in [x + 1, 00) at time t, while J x (t)~~ is the number of particles that lie in [1, 00) at time 
and in (— 00, x] at time t. When more than one process (to, 77, etc) is considered in a 
coupling, the currents of the processes are denoted by J x (t), J x (t), etc. 

Flux and characteristic speed. The average net rate at which particles in the sta- 
tionary (p, g)-ASEP at density p move across a fixed edge (i, i + 1) is the Bux 

(1.1) H(p) = (p-q)p(l-p). 
The characteristic speed at density p is 

(1.2) V» = H'{p) = (p-q)(l-2p). 
In the stationary process the expected current is 

(1.3) EP[J x (t))=tH(p)-xp 

as can be seen by noting that particles that crossed the edge (0, 1) either also crossed (x, x + 
1) and contributed to J x (t) or did not. Another important and less obvious expectation is 

(1.4) W[Q{t)]=tH\p) 

for the second class particle that starts at the origin. We derive (jl.4p in Section [2j 
We can now state the main result, the moment bounds on the second class particle. 

Theorem 1.1. Consider (p, q)-ASEP with rates such that 0<q = l— p<p< 1. 

(Upper bound) For densities p G (0, 1) there exists a constant < C\(p) < 00 such that 
for 1 < m < 3 and t > 1 

(1.5) W{\Q(t)-V p t\ m \ < 

L ' 3 — m 

C\(p) is a continuous function of p. 

(Lower bound) For densities p £ (0,1) there exist constants < to(p),C2(p) < 00 such 
that for t > to(/°) an d all m>l 

(1.6) W{\Q(t)-V p t\ m ) > C 2 (p)t 2m/3 . 
C2(p) and to(p) are continuous functions of p £ (0,1). 

The constants Ci,io>C2 depend also on p but since p is regarded as fixed we omit this 
dependence from the notation. The boundary dependence on p is as follows: as p — > {0, 1}, 
— * 00 ) Cz(p) ~^ and to(p) — > 00. The convergence C^p) — » as p — > {0,1} is 
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natural because for p G {0, 1} Q(t) is a random walk. The divergences of C\(p) and to(p) 
appear to be artifacts of our proof. 

For the third moment and beyond our upper bound argument gives these bounds: for 
< p < 1 

B(p, 3)t 2 log t for m = 3 and t > e, 
B(p, m)t m ~ 1 for 3 < m < oo and t > 1. 



(1.7) E^{ \Q(t) - V p t\ m } < 



The constant B(p,m) is finite and continuous in p for each fixed 3 < m < oo. 
A key identity proved in the next section states that 

w[j VVPn {t)] = p(i- P )w\Q(t)-[vn\\. 

From this and Theorem 11.11 we get the order of the variance of the current as seen by an 
observer traveling at the characteristic speed V p : for large enough t, 

c^p)^ 3 < y a vp[J lVPtl (t)} < c 2 { P )t 2 /\ 

In the case of TASEP (totally asymmetric simple exclusion process with p = 1 = 1 — q) 
P. L. Ferrari and Spohn [5] proved a distributional limit for the current J|y Pi j(i). At the 
time of this writing this precise distributional limit has not yet been proved for ASEP. If 
the observer chooses a speed V ^ V p then Var p [J|y t j (t)] is of order t and in fact J|yt](£) 
satisfies a central limit theorem that records Gaussian fluctuations of the initial particle 
configuration [3]. 

Organization of the paper. Section [2] proves the identities that connect the current 
and the second class particle. In Section [3] we develop the coupling that keeps two second 
class particles in different densities ordered. Section [5] contains a tail bound on a biased 
random walk in an inhomogeneous environment that is needed for the last section. Section 
[5] proves Theorem ll.il 

Further notation. Z + = {0, 1, 2, ... } and N = {1, 2, 3, . . . }. Constants denoted by 
C or Cj (i = 1,2,3, ... ) can change from line to line. Centering of a random variable is 
denoted by X = X - EX. Shift on {0, 1} Z is denoted by (Oi0j)j = uj i+j . 

2. COVARIANCE IDENTITIES 

In this section we can consider more general exclusion processes with bounded jumps 
and rates that depend on the configuration around the jump location. 

Let R be a constant that bounds the range of admissible jumps. For 1 < k < R let and 
qk be functions of particle configurations that depend only on coordinates (coi : —R < i < R) 
and satisfy < pk(oj) , qk{^) < 1- The rule for the evolution is that for all i £ 7L and all 
k G {1, . . . , R}, independently of everything else, the exchange uj \— > u l ' l+k happens at rate 

p k (6iio)uJi(l - io i+k ) + q k (6 i+k Lu)uJi +k {l - uji). 

The configuration u! t,t+k is the result of exchanging the contents of sites i and i + k: 

Wj+fc, ifj = i 
Ui, if j = i + k 
Uj, if j i + k}. 
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The rule of evolution can be restated as follows: whenever possible a particle jumps from i 
to i + k at rate pk{0iU)) and from i to i — k at rate qk(0iuj). 

Key assumption. Bernoulli distributions {^ p }pe[o,i] are invariant for the process. As 
throughout the paper, P p , E p , Var p and Cov p refer to the density-p invariant process. 

We state formulas that tie together the current, space-time covariance (also called "two- 
point function"), and the second class particle. We begin with the well-known formula that 
connects the two-point function with the second class particle. 

Let P p be the probability distribution of two coupled processes uj < co + that start with 
identical Bernoulli- p occupation variables w«(0) = uif(0) at i ^ and a single discrepancy 
Q that starts at the origin. In other words Q(0) = and uj + (t) = uj(t) + Squ\. Then for 
< p < 1 

(2.1) Cov^t), u (0)] = p(l - p)P p {Q(t) = j}. 

Note that the two sides of the identity come from different processes: the left-hand side 
is a covariance in a stationary process, while the right-hand side is in terms of processes 
perturbed at the origin at time 0. For the sake of completeness we give below a proof of 

(EU). 

Theorem 2.1. For any density 0<p<l, zGZ and t > we have these formulas: 

(2.2) Var^[J z (t)] = ^\j - z\ Cov"[w i (t),£«;o(0)] = p(l - p)&\Q(t) - z\. 

and 

(2.3) l-EP[J z (t)]=mQ{t)]-z. 

The terms in the series in (|2.2j) decay exponentially in \j\, uniformly over p. All members 
of these identities are continuous functions of p E [0,1]. At p = and p = 1 the left-hand 
side of (|2.3p is a one-sided derivative. 

Some comments and consequences follow. 

At p = the second class particle Q(t) is a random walk that takes jumps of size k at 
rate Pk($o) and jumps of size — k at rate qk($o)- At the other extreme, p = 1, Q(t) is a 
random walk that takes jumps of size k at rate ^(1 — <5_/%) and jumps of size —k at rate 
Pfe(l — 5k). Here 1 denotes the configuration u = 1. In particular, in (p, g)-ASEP, as p goes 
from to 1, Q interpolates between nearest-neighbor random walks with rates (p, q) and 

(<?,£>)• 

The equilibrium current past the origin satisfies 

(2.4) EP[J (t)} = tH(p) 

where the flux H(p) is the expected rate of particle motion across any fixed edge in the 
stationary density-/) process: 

R 

H{p) = Y J kEP{p k -q k ). 
k=l 
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Combining (12. 3p for z = and (|2.4p gives the useful identity 

(2.5) E p [Q(t)} = tH'(p). 
We turn to the proofs, beginning with (|2.1[) . 

Proof of equation (|2.ip . This is a straight-forward calculation. 

Cov"[c^(t), wo(0)] = E>j(t)wo(0)] - p 2 = ^"[^(t) | wq(0) = 1] - /> 2 
= p[E^[uj{t) | wo(0) = 1] - pEPfait) | wo(0) = 1] - (1 - p)£ p [wj(t) | wo(0) 0] ; 
= p(l - p)(#'[a^(i) | wo(0) = 1] - | £jo(0) = 0]' 

= p(l - p) (p>[w+(t) = 1] - P p h(t) = 1]) = p(l - p)P p [ut(t) = 1 , Uj {t) = 0] 
= p(l-p)P"[Q(t)=j]. □ 

The remainder of this section proves Theorem 12.11 The second equality in (|2.2p comes 
from (gZEJ) . 

Let w be a stationary exclusion process satisfying the assumptions made in this section, 
with i.i.d. Bernoulli(p) distributed occupations {u>i(t)} at any fixed time t. To approxi- 
mate the infinite system with finite systems, for each N 6 N let process uj n have initial 
configuration 

(2.6) uf(0)=Ui(0)l { _ N <i< N} . 

We assume that all these processes are coupled through a Harris-type construction, with 
jump attempts prompted by Poisson clocks, with appropriate rates, attached to directed 
edges for \j — i\ < R. Let J^(t) denote the current in process u> N . 
Let z(0) = 0, z(t) = z, and introduce the counting variables 

(2-7) I»®= £ u£(t), I N (t)= £ ^(t). 

n>z(t) n<z(t) 

Then the current can be expressed as 

and its variance as 

Var jf(i) = Cov(/^ (t) - I* (0), I? (0) - I? (t)) 

= Cov(l£(t),I*(0)) + Cov(J* (()),/*(*)) 

- Cov(/^ (0), 15(0)) - Cov(/^ (i), (*)) 
= £ Cov[^(t),^(0)]+ £ Cov[u,f(t),u^(0)] 

fc<0,m>2 fc<z,m>0 

- Cov(/^ (0), I? (0)) - Cov(/^ (t), (t)). 

Independence gives 

Cov(lf(0),J?(0)) = 
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and the identity above simplifies to 

Varjf(f) = £ Cov[^(t),^(0)] 

k<0,m>z 

2 ' 8 + £ Cov[^(t),^(0)]-Cov(I^(t),/^(t)). 

fc<2, m>0 

We show that identity (|2.8p converges to identity (|2.2p . 

To take advantage of the decaying correlations that result from the bounded jump range, 
we introduce another family of auxiliary processes. Let rj a,b and r] N ' a,b be exclusion processes 
defined from initial conditions 

(2.9) <' fe (0) = Wi(0)l {a<i<6} and 7?f - a ' 6 (0) = u;f (0)l {a<i<6} 

and with a "reduced" Poisson construction: all Poisson jumps that involve any site out- 
side (a, b) are deleted. The point of this definition is that for disjoint intervals (a, b) and 
(u,v), processes r] a,b and rj u,v are independent because they do not share initial occupation 
variables or Poisson clocks. 

We write P for the probability measure under which all these coupled processes {lo,u n , 
r/ a ' b , rj N ' a,b } live. The next lemma is valid for completely general initial occupations {cjj(O)}. 

Lemma 2.2. Let {iOi(0)} be an arbitrary deterministic or random initial configuration, 
and define initial configurations lo n (0), r] a,b (0) and r] N,a,b (0) by (|2.6p and (|2.9p . For a fixed 
< t < oo there is a constant C(t) < oo such that, for all indices N, m, k, and all times 
s G [0, t] 

(2.10) E|w^(s) - Vk ' fe - m ' fc+m (s)| < e~ c ^ m 

(2.11) and EK(s) - ^- m ' fc+m ( S )| < e - c W m . 



Proof. We give the argument for a particular iV and But the argument and the resulting 
bound have no dependence on N, k or the initial configurations. The only issue is the 
disconnectedness of the graph created by the Poisson clocks in the Harris construction 
during time interval [0, t] . 

For any given integer % > 0, there is a fixed positive probability that no site in the 
intervals [k — i — R + 1, k — i] and [k + i,k + i + R — l] is involved in any Poisson jump during 
time interval [0, £]. Hence the probability that this event fails at each integer i = 3R£ for 
1 < £ < n is exponentially small in n. (For distinct £ the clocks involved are independent.) 

If this event succeeds for i = 3R£, then u>j? (s) = 7 ^' fc_m ' fc+m ( s ) f or a ll m > 2>R£ + R 
and s £ [0, t] by the following reasoning. This m is large enough to ensure that the initial 
configurations u> N (0) and ^ Ar : fc_m > fc + m (o) agree on (k — i — R,k + i + R). The absence of 
Poisson jumps has frozen the configurations in (k — i — R, k — i] U [k + i , k + i + R) up to time 
t. Also, up to time t no Poisson jump connects the process inside the interval (k — i,k + i) 
to the outside. The rates of the jumps inside the interval (k — i,k + i) are determined by the 
configuration in the larger interval (k — i — R,k + i + R). Consequently inside (k — i, k + 
up to time t, the processes ui N and ^ JV ' fc -' m ' fc + m execute the same moves. □ 

Now combine the previous lemma with independent initial occupations. 
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Corollary 2.3. Let {w;(0)} be i.i.d. Bernoulli(p) and uj n (0) defined by (|2.6j) . Fix < i < 

oo. There exists a constant C(t) such that for all indices N,i,k, all densities < p < 1 
and times s £ [0, t], 

(2.12) | E5f (s)£f(i) | < 4e~ CW|fc - i| 
and 

(2.13) | Cov"[a;i(a),w fc (t)] | = | B^(s^ k (t) \ < Ae' ^^. 
For any function g on Z £/iaf grows at most polynomially, 

(2.14) [0,1]9^5] 9(k) Cav"[u k (t), wo(0)] 

fe 

zs a continuous function. 

Proof. Let m = |i — &|. Variables r^' 2 m / 2 ,«+m/2^^ 77^^ m/2,fc+m/2^ independent 
and so their covariance is zero. (|2.12p follows from (|2.10p . Similarly for (|2. 13[) . 
For fixed k, m the function 

(2.15) p^Cov"[ V k k - m ' k+m (t),u; (0)] 

is continuous because the expectation depends on p through the finitely many variables 
{r] k ~ m ' k+m (0) : k — m < i < k + m, u>q(0)}. These can be coupled simultaneously for all 
values of p through i.i.d. uniform Ui ~ Unif(0, 1) by writing 

(2.16) ?? fc-m,fc+m,p (()) = 1{Q <Vi< p} 

Continuity in (|2.15p follows by dominated convergence. Estimate (|2.1ip then shows the 
continuity of the individual terms in the series in (|2.14[) . The continuity of the whole series 
follows from the p- uniformity of the tail bounds in (|2.13p . □ 

Once N > \a\ V \b\, r, N ^ h {t) = r] a ' b (t). Fr om this and Lemma 12.21 one concludes that 

(2.17) uof (t) — > u>i(t) as N — > oo, a.s. and in L 2 , for any i. 
From these estimations we derive this lemma: 



Lemma 2.4. The right-hand side of (12. 8p converges as N — > oo to the middle member of 

(E2). 



Proof. We begin by showing the convergence of the last term of (|2.8j) : 
(2.18) Cav(Z?(t),i^(i))->0. 
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This follows from writing, for a fixed K > 0, 

Cov(lf (t), I? (*)) = E E CovWW. (*)) 

fc>2 j<2 

= E E cov(^(t), < 

0<k-z<K -K<j-z<0 

+E E E5f(t)sf(t)+ e E E^wsftf) 

fc>z j<z-K k>z+K -K<j-z<0 

= E E cov^fw, w f 

0<k-z<K -K<j-z<0 

+ 0( E e-°( fc -^+ E e- c (*^>) 

k>z,j<z-K k>z+K,-K<j-z<0 

= E E Cov(^(t), w f(t))+O(e-^). 

0<fc-z</r -K<j-z<0 

For a fixed AT as A" — ► oo 

E E c w (^ (*)> < (*)) -> E E cov( Wfc (t), = o 

0<fc-z<if -K<j-z<0 0<k-z<K -K<j-z<0 

where the vanishing is due to the assumption of i.i.d. variables at time t. Letting K /* cxd 
finishes the proof of (|2.18|) . 

A similar argument, approximating with the variables Tj^' k \ k \/ 2 > k+ \ k \/ 2 ^ j n tails of 
the series and taking the N ^ oo limit in a finite sum, proves 

(2.19) E Cov[£^(t),^(0)]+ E Cov[o;f(t),<(0)] 

k<0,m>z k<z,m>0 

— ► E Cov[cj m (t),w fe (0)] + E Cov[c; fc (t),cj m (0)] 

fc<0,m>2 fc<2,m>0 

= E' n ~~ z\Cov[u n (t),u} (0)]. 
The last equality used shift-invariance. This proves Lemma (|2.4p . □ 



To complete the proof of (|2.2p it only remains to show that the left-hand side of ([27 
converges as Af — » oo to the first member of (|2.2p . The same line of reasoning works again 
because up to an exponentially small probability in m, the current is not altered by removal 
of all jumps that involve sites outside (—m,m). 

Identity (|2.2p is now proved. The statement about the uniform convergence of the series 
and the continuity in p come from Corollary! 



We turn to the proof of (|2.3p . 
Lemma 2.5. For densities < A < p < 1, currents in stationary processes satisfy 

(2.20) E"[J z {t)\ ~ E X [J z (t)} = J" E(i " *) Cov^t), Wo (0)] £». 
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Proof. Let again process with superscript N be the one whose initial configuration is (|2.6p 
with i.i.d. occupations on [-N, N] and vacant sites elsewhere, and P the measure for all 
the coupled processes. Let I N = ^w^(t) be the number of particles in the process u N . 
I N is a Binomial(2A r + l,p) random variable. For < p < 1 

| E [jf (*)] = ± 2 £ ( 2iV m + - P r^-™ E [j N w N 



m 

\ m. I ' 

m=0 
2N+1 

P 1 ~P 



p0--p) 

[recall definitions (|2.7|) ] 



^ P(I N =m){™- 2N +^ m )E[J N (t)\I? 

m=0 

! - jf (t) ■ j>f (o) - p) 



m 



E 



1 Cov[l^(t)-Jf(0), 7^(0) +4(0)] 



P{l-P) 

(2.21) = ( C ov [If (f ) , 7^ (0)] + Cov [j£ (i) - If (0) , If (0)] 

The last equality used Cov[7^(0), 7^(0)] = that comes from the i.i.d. distribution of 
initial occupations. The first covariance on line (|2.2ip write directly as 

Cov[J?(t), I*(0)] = Y, Cov[^(t),^(0)]. 

fc<0, m>z 

The second covariance on line (|2.2ip write as 

Cov[7^(i) - 7^(0) , 7^(0)] = Cov [7^(0) - I N (t) , 7^(0)] 

= -Cov [l N (t), 7^(0)] = - £ Cov[^(t),^(0)]. 

fc<z, m>0 

Inserting these back on line (I2.21h gives 
A E [7f(t)] = — £ Cov[u^(i),u,f(0)] - £ Cov[^(t), «£(())]). 

^ ^ ^ fc<0,m>z fc<z,m>0 

Thus compared to line (12.19P we have the difference instead of the sum. Integrate over 
the density p and take N — > oo as was taken on the line following (I2.19p . This proves the 
proposition. □ 

This lemma together with (|2.ip gives 

(2.22) Ef[J z (t)} - E x [J z (t)} = J\B e [Q(t)} - z) d9 

for < A < p < 1. Couplings show the continuity of these expectations: 

E x [J z (t)] -» E p [J z (t)] and E x [Q(t)] -> E"[Q(t)] as A -> p in [0, 1]. 
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Precisely speaking, if the initial configurations are coupled as indicated in (|2.16|) . the in- 
tegrands converge a.s. to the corresponding integrands. Then bounds in terms of Poisson 
processes give uniform integrability that makes the expectations converge. 

Thus the right-hand side of (|2.22j) can be differentiated in p and identity (|2.3p for < 
p < 1 follows. 

For the one-sided derivatives, consider the case p = 0. If we take A \ in (12. 22ft then 

E p [Jz(t)]= f P (E e [Q(t)]-z)de. 
Jo 

Continuity of the integrand now allows us to differentiate from the right at p = 0. Sim- 
ilar argument for the left derivative at p = 1 completes the proof of identity (12. 3p . The 
continuity of the right-hand side of (12. 3p was also argued along the way. 

3. A COUPLING FOR MICROSCOPIC CONCAVITY 

As observed in (|2.5p the speed of the second class particle in a density-p ASEP is H'(p). 
Thus by the concavity of the flux H a defect travels slower in a denser system (recall that 
we assume p > q throughout). However, in ASEP the basic coupling does not respect 
this, except in the totally asymmetric (p = 1, q = 0) case. To see this, consider two 
pairs of processes (lu + ,uj) and (rj + ,r]) such that both pairs have one discrepancy: u> + (t) = 
u(t) + Squj^ and r] + (t) = rj(t) + 5q v ^. Assume that u){t) > r/(t). In basic coupling the 
jump from state 









"l 


f 











I 








l 







Vi+i_ 











to state 



happens at rate q and results in = i + l>i = Q ri . 

In this section we construct a different coupling that combines the basic coupling with 
auxiliary clocks for second class particles. The idea is to think of a single "special" second 
class particle as performing a random walk on the process of u — rj second class particles. 
This coupling preserves the expected ordering of the special second class particles, hence it 
can be regarded as a form of microscopic concavity. 

This theorem summarizes the outcome. 

Theorem 3.1. Assume given two initial configurations {0(0)} and {0(0)} an d t wo n °t 
necessarily distinct positions (^(O) and Q^(0) on Z. Suppose the coordinatewise ordering 
C(0) > f(0) holds, Q<(0) < Qt(Q), and 0(0) = 0(0) + 1 for i G {Q< (0) , Q* (0)} . Define the 
configuration (~(0) = C(0) - <5 Q c(o)- 

Then there exists a coupling of processes (£~(i), Q^(i),£(t), Q^(t))t>o with initial state 
(£~(0), (0), £(0), Q*(0)) given in the previous paragraph, such that both pairs ((~,Q^) 
and (OQ^) are (p,q)-ASEP's with a second class particle, and Q^(t) < Q^(t) for all t>0. 



To begin the construction, put two exclusion processes £ and £ in basic coupling, obeying 
Poisson clocks {iV^ i=t1 }. They are ordered so that £ > The ( — £ second class particles 
are labeled in increasing order • • • < X m ^i(t) < X m (t) < X m+ \(t) < ■ ■ ■ . We assume there 
is at least one such second class particle, but beyond that we make no assumption about 
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their number. Thus there is some finite or infinite subinterval I C Z of indices such that 
the positions of the Q — £ second class particles are given by {X m (t) : m G I}. 

We introduce two dynamically evolving labels a(t), b(t) £ I in such a manner that X a ^ (t) 
is the position of a second class antiparticle in the ^-process, Xuf\ (t) is the position of a 
second class particle in the ^-process, and the ordering 

(3-1) X a{t) {t) < X b{t) {t) 

is preserved by the dynamics. 

The labels a(t), b(t) are allowed to jump from m to m ± 1 only when particle X m ±i is 
adjacent to X m . The labels do not take jump commands from the Poisson clocks {N l ~* lzil1 } 
that govern (£,C)- Instead, the directed edges (i,i + 1) and (i,i — 1) are given another 
collection of independent Poisson clocks so that the following jump rates are realized. 

(i) If a = b and X a+ \ = X a + 1 then 



(a, b) jumps to 
(ii) If a = b and X a ^\ = X a — 1 then 
(a, b) jumps to < 



(a, 6 + 1) with rate p - 
(a + 1, b + 1) with rate q. 



(a — 1, b) with rate p — q 

(a — 1, b — 1) with rate q. 

(iii) If a ^ b then a and b jump independently with these rates: 



a jumps to 



b jumps to 



a + 1 with rate q if X a+ \ = X a + 1 

a — 1 with rate p if X a ^i = X a — 1; 

6+1 with rate p if X^+i = X^ + 1 

6—1 with rate q if X^i = X^ — 1. 



Let us emphasize that the pair process (£, () is still governed by the old clocks {N l ^ l±l } 
in the basic coupling. The new clocks that realize rules (i)-(iii) are not observed except 
when sites + 1} are both occupied by X-particles and at least one of X a or X^ lies in 

First note that if initially a(0) < 6(0) then jumps (i)-(iii) preserve the inequality a(t) < 
b(t) which gives (|3.ip , (Since the jumps in point (iii) happen independently, there cannot 
be two simultaneous jumps. So it is not possible for a and 6 to cross each other with a 
(a, b) — > (a + 1, b — 1) move.) 

Define processes C~(*) = C(0 — $x a , t \(t) an< l £ + (0 = + ^x bW (t)- I n other words, to 
produce (~ remove particle X a from and to produce £ + add particle X^ to £. The second 
key point is that, even though these new processes are no longer defined by the standard 
graphical construction, distributionwise they are still ASEP's with second class particles. 
We argue this point for (C~,X a ) and leave the argument for (£, X^) to the reader. 

Lemma 3.2. The pair (C~,X a ) is a (p,q)-ASEP with a second class particle. 

Proof. We check that the jump rates for the process (C~ ,X a ), produced by the combined 
effect of the basic coupling with clocks {N l ^ l±1 } and the new clocks, are the same jump 
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rates that result from denning an (ASEP, second class particle) pair in terms of the graphical 
construction as explained in the Introduction. 

To have notation for the possible jumps, let denote an empty site, 1 a ^"-particle, and 
2 particle X a . Consider a fixed pair (i, i + 1) of sites and write xy with x, y S {0, 1, 2, } for 
the contents of sites (i, i + 1) before and after the jump. Then here are the possible moves 
across the edge {i, i + 1}, and the rates that these moves would have in the basic coupling. 



Type 1 


10 - 


-» 01 


with rate p 




01 - 


-» 10 


with rate g 


Type 2 


20 - 


-» 02 


with rate p 




02 - 


-» 20 


with rate q 


Type 3 


12 - 


-» 21 


with rate p 




21 - 


-» 12 


with rate q 



Our task is to check that the construction of (Q~ ,X a ) actually realizes these rates. 

Jumps of types 1 and 2 are prompted by the clocks {iV* — >i±1 } of the graphical construction 
of (£, C)j an d hence have the correct rates listed above. 

Jumps of type 3 occur in two distinct ways. 

(Type 3.1) First there can be a ^-particle next to X a , and then the rates shown above 
are again realized by the clocks {N" 1 ^'^ 1 } because in the basic coupling the ^-particles have 
priority over the X-particles. 

(Type 3.2) The other alternative is that both sites {i, i + 1} are occupied by X-particles 
and one of them is X a . The clocks {N 1 ^ 1 ^ 11 } cannot interchange the X-particles across the 
edge + because in the (£, £)-graphical construction these are lower priority ^-particles 
that do not jump on top of each other. The otherwise missing jumps are now supplied by 
the "new" clocks that govern the jumps described in rules (i)-(iii). 

Combining (i)-(iii) we can read that if X a = % + 1 and X a _i = i, then a jumps to a — 1 
with rate p. This is the first case of type 3 jumps above, corresponding to a ^"-particle 
moving from % to i + 1 with rate p, and the second class particle X a yielding. On the other 
hand, if X a = i and X a+ \ = % + 1 then a jumps to a + 1 with rate q. This is the second case 
in type 3, corresponding to a C^-partide moving from i + 1 to % with rate q and exchanging 
places with the second class particle X a . 

We have verified that the process {Q~ } X a ) operates with the correct rates. 

To argue from the rates to the correct distribution of the process, we can make use of 
the process (C _ , C)- The processes ((~ ,X a ) and £) determine each other uniquely. The 
virtue of £) is that it has a compact state space and only nearest-neighbor jumps with 
bounded rates. Hence by the basic theory of semigroups and generators of particle systems 
as developed in [6], given the initial configuration, the distribution of the process is uniquely 
determined by the action of the generator on local functions. Thus it suffices to check that 
individual jumps have the correct rates across each edge {i, i + 1}. This is exactly what we 
did above in the language of ((~,X a ). □ 
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Similar argument shows that (£, X\,) is a (p, q)-ASEP with a second class particle. To 
prove Theorem O take Q< = X a and = X b . This gives the coupling whose existence is 
claimed in the theorem. 

To conclude, let us observe that the four processes (£, £ + , (~, () are not in basic coupling. 
For example, the jump from state 





Ci+l 




'l 


1" 


c 


Ci+i 







1 




Si+1 




1 







Ci+l. 











to state 



happens at rate q (second case of rule (i)), while in basic coupling this move is impossible. 

4. Biased random walk in an inhomogeneous environment 

This section proves an estimate that is needed for the proof of Theorem 11.11 Let Z(t) 
be a continuous-time nearest-neighbor random walk on state-space S C Z that contains 
Z_ = {..., —2, —1,0}. Initially Z(0) = 0. Z attempts to jump from x to x + 1 with rate 
p for x < — 1, and from s to i - 1 with rate g for x < 0. Assume p > g = 1 — p and let 
9 = p — q. The rates on S \ Z_ need not be specified. 

Whether jumps are permitted or not is determined by a fixed environment expressed in 
terms of {0, l}-valued functions {u(x, t) : x 6 S, < t < oo}. A jump across edge {x — 1, x} 
in either direction is permitted at time t if u(x, t) = 1, otherwise not. In other words, u(x, t) 
is the indicator of the event that edge {x — 1, x} is open at time t. 

Assumption. Assume that for all x & S and T < oo, u(x,t) flips between and 
1 only finitely many times during < t < T. Assume for convenience right-continuity: 

u(x, t+) = u(x, t). 

Lemma 4.1. For all t > and k > 0, 

P{Z(t) < —k} < e~ 2dk . 
This bound holds for any fixed environment {u(x, t)} subject to the assumption above. 

Proof. Let Y(t) be a walk that operates exactly as Z(t) on Z_ but is rectricted to remain 
in Z_ by setting the rate of jumping from to 1 to zero. Give Y(t) geometric initial 
distribution 

P{Y(0) = -j} = TT(j) = (l - ^) for j > 0. 

The initial points satisfy Y(0) < Z(0) a.s. Couple the walks through Poisson clocks so that 
the inequality Y(t) < Z(t) is preserved for all time < t < oo. 

Without the inhomogeneous environment Y(t) would be a stationary, reversible birth 
and death process. We argue that even with the environment the time marginals Y(t) still 
have distribution tt. This suffices for the conclusion, for then 

P{Z(t) < -k} < P{Y(t) < -k} = (q/p) h = exp(Hog < e~ 2dk . 

To justify the claim about Y(t), consider approximating processes Y^ rn \t), men, with 
the same initial value Y^ m ^(0) = Y(0). Y^ m '(t) evolves so that the environments {u(x,t)} 
restrict its motion only on edges {x — l,x} for — m + 1 < x < 0. In other words, for walk 
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y( m )(t) we set u(x,t) = 1 for x < —m and < t < oo. We couple the walks together so 
that Y(t) = Y^ m \t) until the first time one of the walks exits the interval {— m + 1, . . . , 0}. 

Fixing m for a moment, let = sq < s\ < S2 < S3 < • • • be a partition of the time axis 
so that Sj f 00 and the environments {u(x, t) : —m < x < 0} are constant on each interval 
t G [sj,Sj + i). Then on each time interval [sj,Sj+i) Y^ m \t) is a continuous time Markov 
chain with time-homogeneous jump rates 



c(x, X + 1) 

and 

c(x, x — 1) 



pu(x + 1, Si), —m < x < 
p, x < —m — 1 



qu(x, Si), —m + 1 < x < 
o, x < —m. 



One can check that detailed balance ir(x)c(x,x + 1) = ir(x + l)c(x + l,sc) holds for all 
x < —1. Thus 7r is a reversible measure for walk Y( m \t) on each time interval [sj,Sj+i), 
and we conclude that Y^(t) has distribution ir for all < t < 00. 

The coupling ensures that Y^ m \t) — > Y(t) almost surely as m — > 00, and consequently 
also Y(t) has distribution ir for all < t < 00. □ 



5. Moment bounds for the second class particle 
In this section we prove Theorem I l.li 

5.1. Upper bound. We prove the upper bound of Theorem II .11 by proving bounds for tail 
probabilities. We do this first for the right tail of Q(t). Throughout we assume fixed rates 
p > q = 1 — p and abbreviate 

(5.1) 0=p-q. 
Introduce also the notation 

(5.2) V(t) = B p \Q(t) -V p t\. 

Lemma 5.1. Let r > 1 V S\/6. Then for each density < p < 1 there exists a constant 
C(p) £ (0, 00) such that, for t > 1 and u > rt 2 / 3 , 

(5.3) P?{Q(t) > VH + u}< C{p) ( + 4) • 
The constant C(p) is continuous in p £ (0, 1), and lim p \ i QC(p) = 00. 

Proof. Assume for convenience that u is a positive integer. Since [^J > u/2 for u > 1, (|5.3p 
extends from integers u to real u by an adjustment of the constant C. 

Fix a density < p < 1 and an auxiliary density < A < p that will vary in the 
argument. Start with the basic coupling of three exclusion processes uj > lj~ > r] with this 
initial set-up: 

(a) Initially {wi(0) : i 7^ 0} are i.i.d. Bernoulli^) distributed and ojq(0) = 1. 

(b) Initially w -(0) = w(0) - <5 . 
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(c) Initially variables {??i(0) : i ^ 0} are i.i.d. Bernoulli(A) and 770(0) = 0. The coupling 
of the initial occupations is such that u>i(0) > 771(0) for all i 7^ 0. 

Recall that basic coupling meant that these processes obey common Poisson clocks. 

Let Q(t) be the position of the single second class particle between u(t) and u~(t), 
initially at the origin. Let {Xi(t) : i E Z} be the positions of the to — 77 second class 
particles, initially labeled so that 

• • • < X_ 2 (0) < X_i(0) < X o (0) = < Xi(0) < X 2 (0) < 

These second class particles preserve their labels in the dynamics and stay ordered. Thus 
the Lv(t) configuration consists of first class particles (the 77 (t) process) and second class 
particles (the Xj(t) 7 s). P denotes the probability measure under which all these coupled 
processes live. Note that the marginal distribution of (uj,uj~,Q) under P is exactly as it 
would be under P p . 

For xgZ, J£(t) is the net current in the w-process between space-time positions (1/2, 0) 
and (x + 1/2, t). Similarly J%(t) in the 77-process, and Jx~ v {t) is the net current of second 
class particles. Current in the o;-process is a sum of the first class particle current and the 
second class particle current: 

(5.4) J£(t) = J£(t) + •#-"(*)• 

Q(t) is included among the {Xj(t)} for all time because the basic coupling preserves the 
coordinatewise ordering oj~(t) > rj(t). Define the label rriQ(t) by Q(t) = X mq ^(t) with 
initial value 7Bq(0) = 0. We insert a bound on the label. 

Lemma 5.2. For all t > and k > 0, 

P{m Q (t) >k}< e- 29k . 



Proof of Lemma \5.2l In the basic coupling the label mqit) evolves as follows. When X mQ _\ 
is adjacent to X mQ , ijiq jumps down by one at rate p. And when X mQ+ i is adjacent to X mQ , 
itiq jumps up by one at rate q. When X mQ has no X-particle in either neighboring site, 
the label uiq cannot jump. Thus the situation is like that in Lemma 14.11 (with a reversal 
of lattice directions) with environment given by the adjacency of X-particles: u(m, t) = 
l{X m (t) = X m _i(i) + 1}. However, the basic coupling mixes together the evolution of the 
environment and the walk tuq, so the environment is not specified in advance as required 
by Lemma l4.1i 

We can get around this difficulty by imagining an alternative but distributionally equiv- 
alent construction for the joint process (77, u;~, oS). Let (rj,uj) obey basic coupling with the 
given Poisson clocks {N x ^ x±l } attached to directed edges (x,x ± 1). Divide the u> — 77 
particles further into class II consisting of the particles uj~ — 77 and class III that consists 
only of the single particle u — uj~ = 5q. Let class II have priority over class III. Intro- 
duce another independent set of Poisson clocks {N x ^ x±1 }, also attached to directed edges 
(x, x ± 1) of the space Z where particles move. Let clocks {N x ^ x±1 } govern the exchanges 
between classes II and III. In other words, for each edge {x, x + 1} clocks N x ^ x+1 and 
jyx+i-^x are b serve d if sites {x, x + 1} are both occupied by ui — 77 particles. All other 
jumps are prompted by the original clocks. 
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The rates for individual jumps are the same in this alternative construction as in the 
earlier one where all processes were together in basic coupling. Thus the same distribution 
for the process (77, lj~ , u) is created. 

To apply Lemma 14.11 perform the construction in two steps. First construct the process 
(?7,cj) for all time. This determines the environment u(m,t) = l{X m (t) = X m -i(t) + 1}. 
Then run the dynamics of classes II and III in this environment. Now Lemma 14. 1 1 gives the 
bound for ijiq. □ 

The preliminaries are ready and we begin to develop a series of inequalities. Let u and 
k be positive integers. 

P{Q(t) > V p t + u} 

(5-5) < P{m Q (t) >k} + P{J£ VPti+u (t) - J v lVPt]+u (t) > —k}. 

To explain the inequality above, if Q(t) > V p t + u and rriQ(t) < k then AQ.(i) > [V p t\ + u. 
This puts the bound 

Ji \vn\+u^) > ~ k 

on the second class particle current, because at most particles Xi, . . . , X^-i could have 
made a negative contribution to this current. 

Lemma 15.21 takes care of the first probability on line (|5.5h . We work on the second 
probability on line (|5.5p . 

Here is a simple observation that will be used repeatedly. Process u can be coupled 
with a stationary density-p process so that the coupled pair (lo,uj^) has at most 1 
discrepancy. In this coupling 

(5.6) \JZ(t)-jf\t)\<l. 

This way we can use computations for stationary processes at the expense of small errors. 

Recall that V p = H'(p). Let c\ below be a constant that absorbs the errors from using 
means of stationary processes and from ignoring integer parts. It satisfies |ci| < 3. 

E-TfUj+Jt) - EJf v , tJ+u (t) = tH(p) - (H'(p)t + u)p - tff(A) + (H\p)t + «)A + c x 

= -\tH"{p)(p-\) 2 - U (p-X) + Cl 

(5.7) = tO(p - A) 2 - u(p - A) + ci. 

For more general fluxes with nonvanishing H"{p) the Taylor expansion would produce more 
terms above. 

The discussion splits into three cases according to the range of u- values. Only the first 
case requires substantial work. 

Case 1. rt 2/3 <u< pdt. 

Choose 

,2 

3. 



u 

(5.8) A = p — and k 

K ' H 29t [80*. 

The assumptions r > 1 V 8\/#, t > 1 and u > rt 2 / 3 ensure that 



u 
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In the next line below the —3 in k absorbs c\ from line (15, 7h . Recall that X = X — EX 
stands for a centered random variable. We continue bounding the second probability from 
line ([531) . 

< P{^tj + n(*) " Jlvn\ + u® > ^} 

< -^Var{j[V PtJ+u (i) - Jf v „ J+u (t)} 

(5-10) < -^(Var{J^ tJ+u (t)} + Var{ ^ J+B (t)} ) . 

We develop bounds on the variances above, first for J w . Pass to the stationary density-/) 
process via (|5.6p and apply (|2.2p of Theorem 12, It 

Var{j^ J+M (t)}<2Var'{j 

\ypt\+u W} + 2 
= 2p(l - p)E|Q(t) - L^ p iJ - n| + 2 
< E|Q(t) - + n + 3. 

As was already pointed out, as far as Q(t) goes, the E p expectation in the right member of 
()2.2p is the same as E in the coupling of this section. Recall the notation ^(i) from (|5.2I) . 
bound 3 by 3« (recall that u > 1), and write the above bound in the form 

(5.11) Var{jj- VPtJ+u (t)}<*(t)+4u. 

In order to get the same bound for Varj J?yp t \ +U (t)} we utilize the coupling developed in 

Section Let Var A denote variance in the stationary density- A process and let Q v (t) denote 
the position of a second class particle added to a process r\ defined as at the beginning of 
Section 15.11 

Reasoning as was done for (|5. 1 1 [) : switch to a stationary density- A process and apply 
CL2D: 

Var{j^ tJ+u (t)} < 2Var A {J LyPtJ+M (t)} +2 

< E x \Q v (t) - [V p t\ -u\+2 
<E x \Q\t)-V p t\+Au 

Introduce process ((~(t), (t), rj(t), Q v (t))t>o coupled as in Theorem 13.11 where ( starts 
with Bernoulli(p) occupations away from the origin and initially Q<(Q) = Q1(0) = 0. Below 
apply the triangle inequality and use inequality Q<(t) < Q v (t) from Theorem ED Thus 
continuing from above: 

= E|Q"(t) - QHt) + Q C (t) ~V p t\+ Au 

< E{Q"(t) - QHt)} + V\Q C (t) ~V p t\+ Au 
= V x t - V p t + (t) + Au 

= 26t(p - A) + *(t) + Au 
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(5.12) =¥(t)+5u. 

Marginally the process (CQ**) is the same as the process (u>,Q) in the coupling of this 
section, hence the appearance of above. Then we used (|2.5p for the expectations of 
the second class particles and the choice (j5.8[) of A. 
Insert bounds (pUl) and flEIZj) into (IBTTUI) to get 

(5.13) P^Ptj+uW " JfvPtj+uW > ~*> C92 {^ {t) + 
where C = 1152. By ([531) and Lemma 

P{m Q (t) >k}< e-" 2 / 8 * 
which we bound by (recalling u > 1) 

e -« a /M<cr( * y<c4- 

Insert this and (15.13P into line (15. 5p to get the upper tail bound 

2 2 

(5.14) P{Q(t) > F't + u} < C(0) ( V W + t) 
which is exactly the goal (15. 3p for Case 1. 

Case 2. p6t < u < 3t. 

This comes from the previous case. With u < 3t, apply (|5.14p to v = ^pdu < pOt to get 
the bound 

P{Q(i) > V p t + u}< P{Q(t) > V p t + v} 

(5-15) < C(0) ( 4*W + 4) < ( 4*W + 4) • 

Case 3. u > 3i. 

This comes from a large deviation bound. Q(t) is stochastically dominated by a rate 1 
Poisson process Z% and \V P \ < 1. A straightforward exponential Chebyshev argument for 
the Poisson distribution gives the bound 

(5.16) ? p {Q(t) > V p t + u}< P{Z t > 2u/3} < e~ Bu 

for t > 1 and u > 3t with a constant -B independent of all the parameters. The rightmost 
member of (|5.16p is dominated by Ct 2 /u 3 with constant C a fixed number. So it can 
be covered by the rightmost member of (|5.15p by increasing the constant C{6) because 
p~ A > 1. Renaming the constant in (|5.15p gives the conclusion of Lemma l5.11 □ 

We extend the bound to both tails. 

Lemma 5.3. Let r > 1 V 8y/6. Then for each density < p < 1 there exists a constant 
Ci(p) £ (0, oo) such that, for t > 1 and u > rt 2 / 3 , 

2 2 

(5.17) P'{ \Q(t) - V p t\ >u}< dip) ( + -3 ) . 
TTie constant C\(p) is continuous in p G (0, 1), and linip^jo,!} Ci(p) = oo. 
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Proof. Consider coupled processes (u>,u> + ) that start with i.i.d. density p £ (0,1) occupa- 
tions away from the origin and initially w + (0) = w(0) + 5q- So there is one second class 
particle Q initially at the origin. Define new processes 

Qi(t) = 1 - Ui(t) and Q^(t) = 1 - ujf{t). 

Process u records the dynamics of holes in the w-process, and similarly for lD~. The 
discrepancy between (Q~,Q) is exactly the same as the discrepancy between (lo,uj + ). That 
is, Q(t) = Q(t). Processes (Q~,uj) are instances of (p, q)-ASEP with density p = 1 — p and 
rates p = q and q = p. 

To recover the original rates (p, q) we reflect the lattice across the origin. Define 

wftt) = oU(t) and v?-(t) = ul i (t). 

Process (u> R ~,u> ) is an instance of (p, c/)-ASEP at density 1 — p with a discrepancy, so the 
previous bound (|5.15p applies. The discrepancy is now at Q R (t) = —Q(t) = —Q(t). The 
characteristic speed is V x ~ p = —V p . Hence also 

tf(t) = E p \Q(t) - V p t\ = E p \Q R (t) - V l ~ p t\. 

By an application of inequality (|5.15p to the process (ui R ~ , lo r , Q r ), 

P p {Q(t) < V p t -u} = P p {Q R (t) > V l ~ p t + u} 

Combining (|5.15p and (|5.18p we have the conclusion (|5.17p with constant 

(5-19) Ci(p) = a ^ 

V ; VP; p 4 A(l-p) 4 

where C(8) is a constant that depends on 9. □ 

Let us also record the two-sided large deviation bound, one side of which was argued 
above in (I5.16p . 

Lemma 5.4. There exists a constant B such that 

(5.20) P p { \Q(t) - V p t\ >u}< e~ Bu 

for t > 1, u > 3t and < p < 1. 

Inequalities (|5,17p and (|5.20p give the upper bound on the moments of the second class 
particle in (|1.5p via a two-step integration argument. We keep track of the precise constants 
for a while for use in the lower bound proof to come. First by (|5.17p 



POO 

¥(*)=/ P p {\Q(t)-V p t\ >u}du 
Jo 

^*< t >+("+^> 2/s ' 



d.u 
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Fixing r = max{l, 8y/9, Ci(p) 1 / 3 } shows that < C 2 (/o)t 2//3 for a new constant C 2 (p) 
with the same properties as C\{p). 

Put this back into the estimate (|5.17j) to get 

/ t 8 / 3 t 2 \ t 2 

(5.21) P^{ \Q(t) - Vt\ >u}< C 3 (p) ( — + -)< C 3 (p)-3 

\ Uj Hi / Hi 

with new constant C 3 (p). The second inequality used u > rt 2 ! 3 . Now take m > 1 and this 
time use both (foTTTD and (foT2TJD : 

/>oo 

E"iQ(t) - vnp = "» / p p { \Q(t) - v p t\ > u} u m ^ du 

Jo 

r3t poo 

(5.22) <r m t 2m / 3 + C 3 (p) t 2 u m - 4 du + m e~ Bu u m - x du. 

Jrt 2 /3 J3t 

Performing the integrations gives these bounds: 

^(fZL t 2™/3 for 1< m < 3 and t > 1, 
3 — m 

w\Q(t) - vt\ m < i c ^ p)t 2 logt for m = 3 and t > 6j 

CQ(p,m)t m ~ 1 for 3 < m < oo and t > 1. 

One can see that the constants Ci are continuous in p. Cq diverges to oo as m \ 3 or 
m y oo. 

After renaming constants, this gives the upper bounds in Theorem 11.11 and (|1.7f) . We 
record one case here for use in the lower bound proof. 

Proposition 5.5. There exists a continuous function < Cjjb(p) < °° of density < p < 1 
such that for t > 1, 

E'|Q(t) - yt| < C UB {p)t 2 ' 3 . 



5.2. Lower bound. Fix a density < p < 1. For the lower bound we prove that t" 2 / 3 ^(t) 
has a positive lower bound for large enough t where ^(t) was defined in (|5.2p . The lower 
bound (jl.6p follows then for all m > 1 from Jensen's inequality. 
Let ai , a2 be finite positive constants and set 

(5.23) b = o|/(32(9). 

In the argument ai and 02 will be chosen sufficiently large relative to p and relative to the 
constants Cjjb{~) m Proposition 15.51 Define an auxiliary density A = p — bt~ l l 3 . Fix to > 1 
and large enough so that A 6 (p/2,p) for t > to- Restrict all subsequent calculations to 
t > to- Define positive integers 

(5.24) u = Lait 2/3 J and n = [V x t\ - [V p t\ + «. 
Start with a basic coupling of three processes 77 < r/ + < £: 

(a) Initially ry has i.i.d. Bernoulli(A) occupations {^(0) : i ^ —n} and r/_ n (0) = 0. 

(b) Initially i] + (0) = r/(0) + <J_ n . <2 ( ~ n) (t) is the location of the unique discrepancy 
between 77 (t) and r] + (t). 

(c) Initially £ has independent occupation variables, coupled with 77(0) as follows: 



22 



M. BALAZS AND T. SEPPALAINEN 



(c.l) ^(0) = th(0) for -n < i < 0. 

(C.2) C-n(0) = 1. 

(c.3) For i < —n and i > variables Q(0) are i.i.d. Bernoulli^) and Q(0) > rji(0). 

Thus the initial density of £ is piecewise constant: on the segment {— n + 1, . . . , 0} £(0) 
is i.i.d. with density A, at site — n £(0) has density 1, and elsewhere on Z £(0) is i.i.d. with 
density p. 

Label the C ~ V second class particles as {Y m (t) : m £ Z} so that initially 
• • • < K.i(0) < Y (0) = -n = Q { - n) (0) < < Fi(0) < Y 2 (0) < 

Let again toq(£) be the label such that Q( -n )(t) = Y mQ ^(t). Initially mg(0) = 0. The 

inclusion Q(~ n )(i) G {Y m (t)} persists for all time because the basic coupling preserves the 
ordering > r] + (t). Through the basic coupling mg jumps to the left with rate q and to 
the right with rate p, but only when there is a Y-particle adjacent to Y mQ . As in the proof 
of Lemma 15.21 we can apply Lemma [4.11 to prove this statement: 

Lemma 5.6. For all t > and k > 0, 

P{m Q (t) < —k} < e- 29k . 



Fix a\ > large enough so that u = [ait 2//3 J satisfies, for all the auxiliary densities A, 

(5.25) P{Q(~") (t) > [V p t\ } = P{Q(-") (t) > -n + L^ A iJ + u} < |. 

This is possible because the constant C[/b(A) in the upper bound in Proposition 15.51 is 
continuous in the density and we restricted to t > to to ensure that A 6 (p/2,p). 

If Q(~ n \t) < [VH\ and m Q (t) > -k then K_ fc (t) < [V p t\. As was argued for (Q in 
Section I5TT1 this implies a bound on the second class particle current: 

J\ypt\ (*) ~~ ^L^ P *J ^ = "^LV^J — ^" 

So from (fS35|) 

J<p{gt-»)(t)< [y'tj} 

(5.26) , 

< P{m Q (t) < -k} + P{ jf yptJ (t) - J* VPti (t) < k}. 

Take a 2 > large enough so that 



(5.27) a 2 >S + 8y/C UB (a) for a € (p/2, p). 
Increase io further so that for t > to 

(5.28) P{m Q (i) < -a 2 ^ 1/3 + 3} < i 

(utilizing Lemma l5,6j) . Combine displays (|5.26p and (15, 28ft with k = [fl^ J — 2 to get the 
next inequality. Then split the probability. 

\ < ViJfyptl (*) " J [Wj (*) < «2i 1/3 - 2} 
< P{jf v „, (i) < 2a 2 i 1/3 + *0(2pA - A 2 )} 



' \vn\ ' 



(5.29) + P{J* . (t) > a 2 t 1/3 + t0(2pA - A 2 ) + 2}. 
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We treat line (|5.29j) . Recall that P x denotes probabilities of a stationary density- A 
process. As described above (|5.6p we can imagine a basic coupling in which the ^-process 
differs from a stationary density-A process by at most one discrepancy. Compute the mean 
current in the stationary process: 

E x {J [VPtl (t)}=tH(X)-X[V p t\ 

< tH(X) - XV p t + 1 = t9(2pX - A 2 ) + 1. 

From these comes the bound 

line (|fT29j) < P x { J [VPtl (t) > a 2 t 1/3 + t6{2 P X - A 2 ) + 1} 

< P X {J [VPt] (t) > a 2 t 1 ' 3 } < fl 2 2 r 2 / 3 Var A {J L ^j(t)} 

B x \Q(t) - \VH\\ E x \Q(t)-V x t\ 29bt 2 / 3 + 1 
- ^p75 - ^p75 + a 2 t 2 /3 

(5.30) < CuB {X)af + ^ + ^<\. 

After Chebyshev above we used (|2.2p and introduced a second class particle Q(t) in a 
density-A system under the measure P A . [V p t\ was replaced with V x t at the cost of an 
error 1 for dropping integer parts, and V p — V x = 29bt~ 1 / 3 . Last we applied the upper 
bound from Proposition 15.51 and properties (|5.23|) . (|5.27|) and t > to > 1. 
Put this last bound back into line (|5.29p . This leaves 

(5.31) § < P{jf yptJ (t) < 2a 2 t 1 ' 3 + t9(2p\ - A 2 )} 

To treat this probability we take the estimation back to a stationary density-/) process 
by inserting the Radon-Nikodym factor. Let 7 denote the distribution of the initial £(0) 
configuration described by (a)-(c) in the beginning of this section. As before v p is the 
density-p i.i.d. Bernoulli measure. Their Radon-Nikodym derivative is 

/H = ^(u,) = ~i{u,_ n = i}- fl (^i{ Wi = i} + i^i{« i = o}). 

V P i=-n+l P P 

Recalling that n = 0(t 2/3 ) and p - A = 0(t~ 1/3 ), 

(5.32) E p (f 2 ) = -(l + < p -i e -(P-A) 2 /p(i-P) < Cl(p) 2 ) 

p\ p(l-p)/ 

where c\{p) is independent of t, continuous in p but diverges to 00 as p — > {0,1}. 
Let A denote the exclusion process event 

A = {J\yn]{t) < 2a 2 i 1/3 + t6(2pX - A 2 )}. 

Then from (|5.3ip 

§ < P{C G ^4} = J P"{A)~<{du) = J P<"(A)/(u,) v p (dio) 

(5.33) < {P p (A)) 1/2 {E p (f)) 1/2 < Cl (p){P p (A)) 1/2 . 
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In the final calculation we bound P P (A) by Chebyshev to return to the current variance. 
Note that 

EP {j[VPt\(t)} = tH(p) - p[V p t\ = t6p 2 +pV p t - p[V p t\ > t6p 2 . 
Also, if we choose a 2 large enough, namely a 2 > (2O480) 1 / 3 , we can ensure that 
(5.34) c 2 = b 2 9 - 2a 2 > 0. 

From line (|5.33p we have, utilizing A = p — fct" 1 / 3 and (15. 341) . 

(8 Cl (p))- 2 < PP{A) = PP{J\y Hl (t) < 2a 2 t 1 / 3 + t6(2pX - A 2 )} 

< P p {J lVPt ] (t) < 2a 2 i 1 / 3 - t6(p - A) 2 } 
<P p {J\vn\(t) < -{b 2 0-2a 2 )t l l z } 
<c 2 2 t- 2 ^Y^{J lVPtl (t)} 

< c 2 2 t- 2 / 3 V(t) 

by (12. 2p again and the abbreviation (|5.2p . We have assumed t > tQ. The first and last lines 
of the calculation above show that t~ 2 ^^(t) has a positive lower bound for t > to- This 
lower bound depends on p through c\(p) in (15.32h . and vanishes as p — > {0, 1}. 
The lower bound of Theorem 11.11 is proved. 
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